First-principles calculations of the self-trapped exciton in crystalline NaCl 
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The atomic and electronic structure of the lowest triplet state of the off-center (C2v symmetry) self- 
trapped exciton (STE) in crystalline NaCl is calculated using the local-spin-density (LSDA) approx- 
imation. In addition, the Franck-Condon broadening of the luminescence peak and the aig — > h^u 
absorption peak are calculated and compared to experiment. LSDA accurately predicts transition 
energies if the initial and final states are both localized or delocalized, but 1 eV discrepancies with 
experiment occur if one state is localized and the other is delocalized. 
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Unlike a molecule, an extended system such as a solid 
can support both spatially localized and delocalized sin- 
gle particle states and excitations. Physical properties 
such as luminescence can differ dramatically depending 
on which type of state occurs. Deciding theoretically 
whether a localized or delocalized solution exists is a 
challenging problem Here we examine NaCl, a clas- 
sic example [|| where electronic excitations self-localize 
by coupling to the lattice, creating local lattice distor- 
tions. Because the degree of localization will affect the 
Coulomb energies, approaches that incompletely cancel 
the self-interaction contribution to the exchange energy 
(e.g., the local density approximation) sometimes fail to 
predict the actual localized solution. 

The ground state of alkali halides with one electron re- 
moved is the Yk center : the resulting hole does not 
delocalize at the top of the valence band, but rather (sym- 
metrically) attracts two Cl~ ions into a tightly bound 
molecule [HH], effectively becoming a Cl2^ molecular ion. 
The local symmetry of this atomic configuration is D2h- 
An excess electron in bulk NaCl forms a large mobile 
(Frohlich) polaron but in the presence of a self- 
trapped hole forms a self-trapped exciton (STE). It has 
been suggested that the self-trapped exciton state 
breaks symmetry and sits off-center with C2v symmetry. 

Although there have been many theoretical studies 
of the self-trapped exciton and Vk center prob- 
lems, no density-functional calculations have been re- 
ported. In this letter we report local-spin-density approx- 
imation (LSDA) calculations of the STE and Vk center. 
The advantage of LSDA calculations is that they provide 
one of the simplest tools capable of providing a realis- 
tic model of this competition. In both cases, our LSDA 
calculations give undistorted, delocalized solutions with 
lower energy than the self-trapped solution, contrary to 
experiment. For the Yk center, no metastable local min- 
imum trapped solution was found; however, for the (neu- 
tral) STE we find locally metastable solutions, with the 



on-center STE 0.14 cV higher in energy than the off- 
center STE solution, which in turn is higher by 0.20 eV 
than a free electron-hole pair. Even with this discrepancy 
in the total energy, the atomic positions for the STE so- 
lution (see Table |) are reasonable and agree well with 
Hartree-Fock second-order M0ller-Plesset (MP2) pertur- 
bation theory ||l^ . We focus on the properties of the local 
minimum solution for the off-center STE, which provides 
a test of the ability of density functional methods to treat 
the strong coupling of electronic and lattice degrees of 
freedom, and calculate the spectral properties of excited 
states. 
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FIG. 1. The dotted contours show the total valence charge 
density in the STE. The solid contours show I^Pp for the 
trapped (a) hole and (b) electron. The open (filled) squares 
represent displaced (ideal) positions of the Na atoms; filled 
circles represent undisplaced CI atoms. The lowest contours 
are 0.04 for the total charge, 0.01 (0.002) for the hole (elec- 
tron) state, with increments of 0.05, 0.01 and 0.002 e (a.u.)~^, 
respectively. 

To solve the LSDA equations we use a plane wave 
pseudopotential method |12 14| with a spin-dependent 
exchange-correlation potential |p^ , and full structural 
relaxation in a supercell approach. For most calcula- 
tions, we used a 32 atom supercell with translation vec- 
tors (2,0,0), (0,2,0) and (1,1,1), giving a nearest neighbor 
distance between STEs of 9.4 A, and used four special 
k-points in the irreducible wedge for the Brillouin zone 
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integrations. Tests varying the number of k-points and 
supercell size suggest that these parameters are adequate 
pGf . One of the in x — y plane nearest CI and two Na 
atoms are most displaced while the rest of the atoms 
move by a much smaller amount. Some calculated struc- 
tural parameters for the STE are given in Table |. 

The calculated atomic displacements from the ideal 
NaCl structure for the off-center STE are shown in Fig. 
|l|. The bsu hole state in the triplet STE is localized on 
the Cl^ molecule (Fig. |l|a), with nearly equal weight on 
the two CI ions. The last spin up electron (aig) is mostly 
localized on the (1/2,1/2,0) vacant halogen site (Fig. |l|b), 
as in the case of the F-center. The formation of the Clj^ 
"molecule" in the STE is mainly due to the shift of a 
single CI. This asymmetric shift can be rationalized by 
noting that the Madelung energy (with canonical charges 
of ±1 for Na and CI) of the D2h configuration is about 
0.17 eV higher than the C2v one, i.e., the ionic Madelung 
terms that favor the rock salt structure in the first place 
favor keeping one of the CI ions on a lattice site; in addi- 
tion, the extra electron in the STE (compared to the Vk 
center) can lower its energy by this distortion. 
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FIG. 2. (a) density of states (states/eV spin formula unit) 
of pure NaCl crystal, where the two spin orientations are 
equivalent. IPR of triplet STE for (b) spin up and (c) spin 
down electrons. The state labeled ba^ is empty and the state 
aig is occupied in the STE. The D2h representations hsu and 
aig became indistinguishable (both Ai) in the true C2v sym- 
metry of the off-center state. 

To obtain vibrational properties of the off-center STE, 
we made finite displacements from the equilibrium geom- 
etry. The Cl^ stretching mode Wstr = 242 cm~^ is smaller 
than the experimental Raman frequency 361 cm^^ 
/,From the force matrix associated with the Cl^ stretch- 
ing mode, we found that only the two neighboring Na 
ions at (5, 1,0) and (1,^ ,0) couple significantly. Unhkc H 
center calculations ||l^, coupling to the Na atoms yields 
only a small shift in the frequency of the Cl^ mode to 
LUstr = 234 cm" ^ Since the LSDA places the STE too 
high in energy, it is not too surprising that the curvature 



of the STE local minimum is underestimated. 

Figure || shows the computed density of states (DOS) 
of the perfect NaCl and the inverse participation ratio 
(IPR=aQ/8 / |^'i(r)|^c?'^r) for states of majority and mi- 
nority spin for the off-center STE. (The IPR is a measure 
of the localization of a state.) The localized hole ba^ and 
electron aig states lie in a gap of about 6 eV between 
the conduction and valence bands of the perfect crystal, 
where as usual the LSDA underestimates the gap. Rather 
than using the single-particle eigenvalues, we obtain esti- 
mates of the excitation energies as the difference between 
total energies of different electronic configurations. To 
calculate the energy to create a free electron-hole pair 
(the gap energy), we occupy spin up states with one ex- 
tra electron, while spin down states have one empty state. 
The energy difference between the two solutions (6.44 eV) 
for the same atomic positions should correspond to the 
free electron- hole pair (experimental value 7.96 eV |p^). 
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FIG. 3. Potential curves for pure NaCI (a), self-trapped 
exciton in the ground state (b) and electronically excited state 
(aig — > bsu) (c). The coordinate a measures the magnitude 
of the displacement, such that q = and a — 1 correspond 
to the undisplaced NaCI crystal and the relaxed STE atomic 
configurations, respectively. 

In the distorted STE solution, the energy of electron- 
hole pair recombination was found by comparing ener- 
gies of the STE solution with the energy of NaCl hav- 
ing the same atomic displacements as that of the STE. 
This energy of 4.25 eV is roughly the same as the en- 
ergy difference between aig and h^u states given by the 
DOS, and compares to a value of 3.35 eV obtained from 
a luminescence experiment [ pO[ at 11 K. The distorted 
NaCl, with a lattice distortion energy of 2.5 eV relative 
to the ideal NaCl positions, is in a highly excited vi- 
brational state. Thus, significant Franck-Condon effects 
in the spectral properties of the STE are expected. To 
qualitatively describe the luminescence, a ground state 
potential curve was calculated for the configuration co- 
ordinate a. It was assumed that all atoms move back to 
the perfect crystal positions proportionally to their dis- 
tortions in the STE solution. The result is shown in Fig. 
^, along with a quadratic fit kaO? /2 with one adjustable 



2 



parameter ka ■ The effective one-dimensional Schrodinger 
equation for '^(a) describes a harmonic oscillator with 
frequency uta — \/Khli = 109 cm^^. The moment of in- 
ertia / = X] MnSR^ was chosen so that 10^/2 equals the 
kinetic energy of the atoms with mass M„, when they 
move from the initial displacements 5Rn- 

The same type of the potential energy curve was cal- 
culated for the STE. The quadratic fit works only in the 
close vicinity of the exciton metastable minimum. The 
resulting frequency is Wf, = 123 cm~^. The experimen- 
tal temperature T=ll K justifies a zero temperature 
approximation (the STE initial state is the vibrational 
ground state). Since the luminescence peak position cor- 
responds to the vibrational level ?t,«170 of the electronic 
ground state, quasiclassical wavefunctions were used in 
the numerical integral evaluation. The sequence of vibra- 
tional sidebands should be replaced by a sequence of con- 
volved densities of phonon states D{uj). We approximate 
this by a Gaussian, D{uj) exp(— a;^/27^)/-\/27r7 with 
the width 7=43 cm~^ chosen such that the first three mo- 
ments coincide with the experimental pl| phonon DOS. 
This gives a luminescence width of 0.43 eV, while the 
experimental width is 0.63 eV. 
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4. Optical conductivity of self-trapped exciton in 
NaCl normalized per volume of a single exciton ao/2. (a) 
oiiuj): El. to the molecular axis in the x — y plane; cr2(tj): E 
along the z direction; and 0-3 (cij): E along the Cl^ molecule, 
(b): average conductivity compared with experiment p2) . 

The optical response <j{uo) of the long-lived triplet STE 
also has been measured [^,^ . The diagonal part of the 
optical conductivity tensor is: 



ne^N 



E,k) 



(1) 



where /^g is the occupancy of the state \£k). The spin 
state index is included in the band index £. Integration 
over the zone has been performed using 26 k-points in the 



irreducible zone. In Fig. absolute optical conductivity 
curves are shown for the three polarizations, E || (1,-1,0), 
(0,0,1), and (1,1,0) (parallel to the the Cl^ molecular 
axis). Fig. ^ shows the average conductivity a{uj) = 
Xj tTi(w)/3, which is compared with experiment p2| , p3[ , 
rescaled so that the total weights under the both curves 
are the same. 

The first peak in Fig. ^ is centered at 0.95 eV; the 
splitting between the peaks for the three different po- 
larizations is not resolved in the calculations. Most of 
the weight in these peaks comes from transitions of the 
last localized spin-up aig electron into empty conduction 
band delocalized states (see Fig. ^d). 

The second peak centered at 3.58 eV is the aig — > b3„ 
transition for the spin down electron (see Fig. ||c). The 
energy difference between the ground state and electron- 
ically excited exciton state with the same atomic config- 
uration turns out to be the same as the eigenvalue dif- 
ference of bsti and aig states of the ground STE. The ex- 
cited exciton will lower its energy by moving atoms back 
to the undistorted positions of perfect NaCl. To apply 
the Franck-Condon principle, we repeat the same type of 
calculations for the excited exciton as we did for NaCl 
and the ground state STE (Fig. ^). When two CI atoms 
move away from each other, the aig empty state merges 
with the valence CI 3p band, which makes it very difh- 
cult to choose which state to depopulate during the iter- 
ations. Instead we used the results for the ground state 
STE to obtain the energy of the electronically excited 
state by adding eigenvalue difference A^i — X12 between 
two states for which dipole matrix element | {£l\p3\£2) p 
is the largest. Results are shown on Fig. |^ along with a 
quadratic fit (ti;c=116 cm~^) which works for the entire 
range of the parameter a. The delta-function of Eq. (|l|) 
corresponding to the aig^bsu transition was replaced by 
a sequence of convolved Gaussian peaks. 

The transition energy (3.58 eV) between the two lo- 
calized states aig^bsi, (see Fig. ^c) agrees well with the 
experimental peak at 3.8 eV. But when initial and final 
states have different degrees of localization, errors of an 
eV in the luminescence (LSDA: 4.25 eV, expt. H]: 3.35 
eV) and in the optical excitation of the bound electron 
into the conduction band (LSDA: 0.95 eV, expt. [p3| : 
1.95, 2.13 and 2.00 eV for three different field polariza- 
tions) occur. A possible explanation for this error is that 
in LSDA the incomplete cancellation of the large repul- 
sive self-interaction energy J drdr'ni{'r)ni{r')/\f — r'\ of 
a localized state rii — I'I'ip and the corresponding ex- 
change term may destabilize a localized solution in favor 
of a delocalized one. Corrections such as self-interaction 
corrections (SIC) or LDA+U |^ may reduce the er- 
ror; in SIC, shifts in the energy of a state i on the order 
of g are expected. If the degrees of local- 
ization of the two states are similar, however, then the 
LSDA transition energies are reasonable. 
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In summary, we have presented LSDA calculations for 
the STE in NaCl, including the coupling between the lat- 
tice and electronic states. The off-center STE is found 
to be more stable than the on-center STE, but both are 
metastable compared to free electron-hole pairs. Both 
luminescence and optical conductivity, including vibra- 
tional Franck-Condon effects, were also calculated. The 
hole state in the off-center STE is found to be rather 
evenly split between the two CI atoms, but the electron 
state is localized to the vacant site left by the shifted CI 
ion. The density functional description of electronic tran- 
sitions between localized states, such as the aig h^u 
absorption peak, agree well with experiment. For tran- 
sitions between the localized and delocalized states, dis- 
crepancies of order 1 eV with experiment arise. Although 
the LSDA can capture many features of the STE states, 
when a localized solution competes with a delocalized so- 
lution, the incomplete cancellation of the self-interaction 
may destabilize the localized solution. Given the useful- 
ness of the LSDA method for unraveling complex mate- 
rials, it is important to test and develop approaches that 
can treat localized and delocalized states on the same 
footing. 
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TABLE I. The equilibrium Cl-Cl distance rci-ci, displace- 
ments of the nearest Na ions Ajva in (A), for the STE. Present 
results are compared to the Hartree-Fock (MP2) theory |ic[ ]. 



Configuration 


rci-ci (MP2) 


Ajv„(MP2) 


STE Dah 


2.73 (2.654) 


0.35 (0.383) 


STE Czv 


2.59 (2.525) 


0.51 (0.488) 
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